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12.  ABSnUCT  (Uutfi,:m}OOwcnk) 


Cyclostationary  processes  are  an  important  class  of  nonstationary  pro  esses.  In  this  report  we 
consider  nonparametric  estimation  of  the  cyclic  cross-spr.ctrum.  A  periodogram  jased  estimator  is 
studied  and  its  asymptotic  behavior  characterized.  This  extends  the  recent  univariate  work  of 
Dandawate  and  Giannakis  to  the  multivariate  case.  T  ?  results  are  useful  for  a  variety  of  multi-sensor 
cyclostationary  signal  processing  scenarios,  such  as  bearing  estimation. 
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1.  Introdut  tion 


In  this  report  we  study  the  statistics  of  a  periodogram  based  estimate  of 
the  cyclic  cross-spectrum.  Cur  multivariate  work  is  based  heavily  upon  the 
univariate  work  of  Dandawate  and  Giannakis  [1,2],  in  which  they  devel- 
ojjed  the  theory  of  periodogram- based  estimates  of  the  cyclic  spectrum  of 
a  single  cyclostationary  process.  Their  work  extended  conventional  spec¬ 
trum  estimation  theory'  for  stationary  time  series  to  the  cyclostationary 
case,  especially  rely  ng  on  ihe  exposition  of  Brllinger  [3]  and  Brillinger 
and  Rosenblatt  [4,51.  simply  extend  the  work  of  Dandawate  and 
Giannakis  to  the  multivariate  case.  The  present  work  can  thus  be  viewed 
as  c  generalization  of  the  muitivanate  cross-periodogram  spectrum  estima¬ 
tion  theory  for  stationary  processes  as  put  forward  by  Brillinger  [3], 

The  pi  :mary  motivation  for  extending  to  the  multivariate  case  is  to  handle 
multiple  time  series  problems,  as  might  anse  in  multi-sensor  signal  proc¬ 
essing  scenario.-;  Because  of  the  dimensionality  difficulties  involved  we 
focus  on  the  second  order  cross -spectrum  case  and  do  not  consider  esti¬ 
mates  of  higher  order  cyclic  cross  cumulants  or  polyspectra,  while  noting 
that  the  univanate  cyclostationary  theory  has  been  extended  to  include 
arbitrary  Aah  order  i  yclic  s?.ui‘,(ics  [  1 ) 

Cyclic  (cross)- spectrum  estimates  based  on  the  periodogram  are  particu¬ 
larly  appealing  due  to  their  nonparametne  nature,  and  Ixjcau^c  they  allow 
the  use  of  the  fast  Fourier  liansform  (FFT)  algorithm,  which  speeds 
computation  It  is  shown  that  such  estimetes  of  the  i  yclic  cross-spectrum 
art  consistent  and  psympioticaJly  normally  distributed  under  mild  condi¬ 
tions  on  the  rime  series.  This  will  allow  for  optimal  criteria  to  be  devel¬ 
oped  for  detection  and  estimation  schemes  based  on  the  cyclic  cross- 
spectrum.  For  example,  the  optimal  estimation  of  the  time  difference  of 
arrival  between  two  scesors  is  key  for  bearing  estimation  schemes,  and 
this  delay  can  be  estimated  for  cyclostationarv  signals  via  the  cyclic  cross- 
specaaim.  Thus,  the  results  cT  this  report  will  prove  useful  in  the  analysis 
of  such  multi- sen.;or  sccnanos. 


2.  Background 

Consider  the  vector-valued  time  scries  X(0  for  t  =  0,  1,2 . with  real¬ 
valued  components  x,(f),  t  =  1 . These  might  arise  as  the  outputs  of 

multiple  sensors,  for  example.  We  assume  without  loss  of  generality  that 
£lx,(f  )1  =  0  for  each  choice  of  i.  With  x^{t)  and  x^it )  components  of  X(t), 
we  define  their  time-varying  cross -covariance  to  be 

r,_,(t;r)  =  Zr[j,(t  -f-  r)r)(01- 
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If  x^it)  and  x^{t)  are  stationary,  then  r)  =  Cj^fr)  is  the  conventional 
cross-covariance  and  exhibits  no  time  dependence.  Backgi  ound  on  the 
stationary  cross  covariance  can  be  found,  for  example,  in  Priestly 
[6,  chapt..  9],  If  x,(0  and  Xj(t)  are  cyclostationary,  then  r)  is  periodic 
(or  almost  periodic)  in  i,  so  for  each  fixed  lag  r  the  time-varying  covari¬ 
ance  can  be  expanded  in  a  Fourier  series  (e.g.,  see  Besicovitch  [7]).  The 
resulting  Fourier  scries  coefficients  are  the  cyclic  cross-covariance,  given  by 


A 


C,(a:r)=  ! 


ini 

T  — »oo 


1 

t=Q 


(2) 


with  ft  called  the  cycle  frequency.  In  analogy  with  the  stationary  case,  we 
may  also  consider  the  Fourier  transform  of  r)  and  C,j(ft;  r)  with 
respect  to  r,  yielding,  respectively,  the  time-varyirif’  spectrum 

•f=— DC  ' 


and  the  rycUc  spectrum,  given  by 

rrr  -nc- 

To  obtain  the  "aulo”  versions  of  these  functions  wc  simply  set  i  =  j  so  that, 
for  example,  c„(ri  r)  f,(<;  r)  is  the  time-varying  autocorrelation. 

Note  that  C,j(0;  r)  is  the  correlation  as  conventionally  defined  for  station¬ 
ary  processes.  Similarly,  5,_j(0;  a;)  is  the  power  spectrum,  defined  as  the 
Fourier  transform  of  f ’,^(0;  r).  Thus,  the  statistics  of  stationary  processes 
arc  a  special  case  of  cyclostation  ary  statistics  with  cycle  frequency  a  =  0, 
which  corresponds  to  periodic  statistics  with  penod  zero  Note  also  that 
t)  and  S^fa.  ,j)  arc  iimc-invariant  quantiucs,  which  is  a  fundamen¬ 
tal  reason  why  cyclic  statistics  are  an  impsirtant  tool  for  analysis  of 
nonstationary  processes. 

In  this  report  we  analyze  the  propenies  of  a  periodogram-ba  ed  estimate 
of  S,^{a,  ,c).  We  make  use  of  the  finite  Fourier  transjorm  of  x,(0.  defined 
as 


I  1 
y  x,i 
(=0 


and  the  cyclic  cross  perio(L>i(ram,  defined  as 


.iT),  ,  1  v-i7  >  ,  , 

(ft:..  )  --  -  A,  (^  :’(A^ 


(6i 


where  *  denotes  complex  conjugate.  Since  .r,(t)  is  real  valued,  this  may  be 


wntten 


llpia]Lu)  =  ’(U’)x^'>(a  -  cv).  (7) 

The  definition  of  equation  (6)  is  a  generalization  of  the  un  v^ariate  version 
of  the  definition  used  by  Hurd  [8]. 

In  deriving  the  covariance  of  uj),  we  make  use  of  higher  order 

cumulants.  All  of  the  above  definitions  may  be  generalized  to  higher  order 
cumulant  cases.  The  “auto"  higher  order  versions  are  defined  by 
Dandawate  ( 1].  Under  the  assumption  that  <  oo,  we  define  the 

joint  cumulant  of  order  k  of  X(0  as,^ 

.  ,a,(U,-  •  •  ,4) 

for  A:  >  2.  In  the  case  of  stationary  processes 

t’lii .  ,ai,  (I  l  ;)■ 

while  in  the  time- varying  case,  in  analogy  with  equation  (1),  we  explicitly 
maintain  the  time  dependence  as 

Ca,  ai,{l  I  ■  '  '  '  •  I  k)  —  a\ .  ..i,  (  f  ^  I  •  '  •  1  )  ■  (10) 

For  X(t)  cyclostationary'  then,  in  extension  of  equation  (2),  we  define. 


'(7’)/,  AviT). 


.av(«;T,,  -.rt_,)  =  lini  Tl.-  ■ 

I  7“ 


to  be  the  klh-order  cyclic  cumulant  of  X(U  at  cycle  frequency  a.  The 
time-varying  cumulant  spectrum  (poly spectrum)  is  wiitten 


( t;  uj\ ,  ■  ■  ■  .  uJk- 


XI  f'oi.  ■  ■  ■  .Ti-  i)c 


r,  =  -  c  TJ  ,  =  -  -x 


and  the  cyclic  cumulant  spectrum  is  defined  as 


‘“^ai .  .ait  1  ^  1 1  '  ‘  '  ^k-  \} 


.  1 )( 


-  J(u.'|  TJ  +  -  1  T*  _  I  1 


*Scc  Brillingcr  [3,  p.  21 ).  We  use  the  notation  )  to  denote  arbitrary  higher  order  cumulants  of  *ah  order, 

rather  than  using  (•).  For  fixed  order  we  retain  the  previous  notation,  so.  for  cxiniplc,  wc  use  c,j{-)  rather  than 
Ca  ^  )  in  the  socondxKdcr  case. 


3.  Cyclic  Cross-Periodogram 

In  this  section  we  explore  the  cyclic  cross-periodogram  u), 

defined  by  equation  (6),  as  a  potential  estimate  of  tlte  cyclic  crcss- 
spectrum  5, .(a;  uj).  The  asymptotic  propenies  of  [a;  oj)  are  suinma- 
rized  in  a  theorem,  which  shows  that  lj)  does  not  result  in  consis¬ 

tent  estimates  of  S,^(a;  uj).  This  is  analogous  to  the  stationary  case,  where 
the  conventional  cross-periodogram  results  in  inconsistent  estimates  of  the 
power  spectrum,  e.g.,  see  Priestly  [6,  chapt.  9].  As  is  well  known  in  the 
stationary  case,  the  key  to  obtaining  consistent  estimates  is  smoothing  the 
periodogram.  As  shown  by  Dandawate  and  Giannakis,  this  idea  extends  to 
the  cyclic  case  as  well  1 1,2]. 

To  begin  our  study  of  a,’),  we  st  de  assumptions  on  the  time  series 

X(/V  These  mixing  conditions  insure  that  samples  of  X(0  that  are  well 
sepaiated  in  time  are  becoming  statistically  independent.  For  a  general 
discussion  of  mixing  conditions,  see  Brillinger  [.?,  chapt.  1). 

Assumption  la 

Ir.lICa,,.  <  oc.  for  t  =  1 .  ■  ,  A*  -  l.Vt.Vfc. 


Assumption  lb 


Y,  l'^llG;(f:  ■^)i  <  oc.  Vf 

=  -OC 


Assumption  2a 

-.rk:'  <  x .  for  ;  -  1 .  ■  ■  ■  ,  A:  -  1 ,  Vq,  VA:. 

'"1  Tk  1=-SX- 


(16) 


Assumption  2b 


<  X. 


.'O. 


Note  that  Assumptions  lb  and  2b  arc  the  second-order  cases  (A:  -  2)  of  the 
moie  general  Assumptions  la  and  2a,  respectively.  For  clarity  we  have 
included  Assumptions  lb  and  2b  .so  wc  can  refe'^  specif'cally  to  the  sec¬ 
ond-order  case  when  desired.  Wc  proceed  by  stating  a  lemma  that  charac¬ 
terizes  the  rate  of  convergence  tas  T  —*  oo)  of  the  cyclic  cumulant  spec¬ 
trum.  In  the  following  the  notation  =  0(6„)  means  that  la„/6„l  is 
bounded  for  n  sufficiently  large,  e.g.,  see  Brillinger  [3,  chapt.  ?|. 
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Lemma  1  If  Assiunjnion  la  holds,  thcti 


(T-n 

’^ai,-,aij(L  •  ■  •  1  1  )  ^  ^  (-  a  I  (L  '  1  '  '  j  ‘^A;  —  1  )  ^ 

r|=--(7  11 


tr-ii 


- I'^Ar  0| 


=  0{T-'). 


□ 


(18) 


For  tlie  case  of  k  -  2.  Lemma  1  reduces  to  the  following  corollary. 
Corollary  1  If  Assiwipiion  lb  holds,  then 


7=:-(7-1) 


□ 


(19) 


The  next  lemma  describes  the  rate  cf  convergence  of  the  finite  Fourier 
transfomi  of  the  Atth-oider  time-varying  eumulant  spectnim.  We  are  gener¬ 
alizing  Theorem  4.3.2  of  Brillinger  [.3].  The  pi  oof  is  similar  to  that  used  by 
Dandawate,  which  is  a  generalization  of  the  stationary  case  of  Brillinger 
and  Rosenblatt  [4]. 


Lemma  2  If  assumption  2a  holds,  then 


T-l 


•  • .  ,  ^  5a,,..,a,(C;  a;,,  . . . ,  +  0(1). 


t=Q 


□ 


(20) 


For  the  case  of  fc  =  2,  lemma  2  reduces  to  the  following  corollary. 
Corollary  2  If  assumption  2b  holds  then, 


catnixPiuto),  Xp{u;,)}  ^  V  -b  C2(l).  □ 

t=o  ’ 

We  are  new'  prepared  to  state  the  main  result  of  this  section  in  the  follow¬ 
ing  theorem  that  describes  the  asymptotic  bias  and  variance  of  die  cyclic 
cToss-periodogram  3’he  proof  of  tlie  theorem  relies  on  the 

above  lemmas. 


Theorem  1  (Bias  and  \  riance  of  the  cyclic  cross-periodogram). 

Suppose  X(0  is  cyclostationary  with  zero  mean  as  defined  in  section  2.  If 
tO],  Uj,  wy)  exists  and  is  finite,  and  Assumptions  1  and  2  hold,  then 

S,, (a, E\iP  (a;  u;)] ,  ^22 


9 


and 

Vun  ti)]  /i,^!)S,t(a  u.-’ t  //  d\i\  l/.’) 

T— ^  ' 

I  4  /'  . -  //;a  ^’)-  □ 

It  is  apparent  from  equation  (22)  that  u.')  is  asymptotically  un 

biased.  However,  from  equation  (24)  we  conclude  that,  in  general, 

u;)  has  non-zero  covariance  and  is  therefore  an  inconsistent  esti¬ 
mate  of  5>,_j(o  w).  A  consistent  estimator  will  subsequently  be  obtained 
based  on  a  smoothed  version  of  '((y;  a.’).  Before  we  do  ibis  we  consider 
some  special  cases  of  equation  (23)  in  the  next  section,  and  m  particular 
we  develop  the  asymptotic  vaiiancc-covariance  matrix  for 


4.  Special  Cases  and  the  Variance-Covariance  Matrix 


In  this  section  we  consider  some  special  cases  of  the  results  of  Theorem  1. 
First,  we  give  two  useful  symmetry  properties.  As  we  have  shown  in 
theorem  1, 

Sij  (a;  U-)  =  (q;  a>)] ,  ^24 

sc  that 


Sij(a;u>)  =  -  cu)]. 


(25) 


From  equation  (25),  and  using  the  fact  that  the  x,(t)  are  real-valued,  we 
can  deduce  the  following  symmetry  properties: 

Sij(Q;a-u,’)  =  Sj,(a;u),  (26) 

S;,(a;u)  =  (27) 

Next  we  consider  some  important  cases  of  equation  (23).  By  setting  i~j- 
k  =  I  and  using  equation  (26)  in  equation  (23),  we  obtain  the  univaiiate 
results  of  Dandawate  and  Giannakis  (compare  with  Dandawate  [1,2]), 

lim  cov{l'P{a;  u),  I-PiP',  /r)}  =  5,i(u;  -  n,  -n)Sii{Q  -  cu  4-  /x  -  /?;  /r  -  0) 

I  -»oo 

4-  Sii{uj  + fi-~  0)Su{a- u  -  (28) 


Another  important  case  that  may  lie  obtained  from  equation  (23)  is  the 
covariance  of  the  conventional  cross-periodogram  for  stationary  processes. 
From  equation  (23),  we  have 


iiin 

r-.oc 


cot'[pP{0\i 


(29) 
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Note  that,  tor  stationary  pr(x:esses,  l^')  t  0  for  o  ^  0.  Ip  this  ease 

equation  (2^))  reduces  to 


lint  cor{!^P{0-j),li:'{0:u:)} 


('1, 


I  <4.(0; a'),  u- /  0,  ,±7r 

I  5,4(0;u')-h’j,(0;a')  +  u;  =■  0,  ±7r. 


This  expression  i  an  be  found,  for  example,  in  IMestly  [6,  chapt.  9j. 

For  periodically  correlated  (PC)  pnxesses  we  ate  restricted  to  the  set  of 

cycle  frequencies  given  by  o  g  {27rA,77o},  for  k  -  0.  1 . Tq-  1.  Using 

(2.^).  this  leads  to 


1 


I 


lirn 


m 


,271  A:  . 

( 'Tfr-  i  .  fy  ( ~7j^  ■  }  - 

-0  1 0 


^7t  A: 


.s',4(i'.^')5ji(0;(>  -  u), 

^,4(0;  n  -  U)  +  -  u-), 


,  _Z 

^  7h 


rt-* 

7j' 


(31) 


for  n  integer. 

To  complete  our  results  on  the  covariance  of  7,^^- (a;  u;).  we  derive  us 
variance-covariance  nialrix.  From  t  ]uation  (2.3)  it  follows  easily  that 


lim  cov{rl[\a\ijj),I^ki^ 


(a;ij)}  =  5i4(();u7)5’j/(0;a  -  u))  -I-  Su{2u 


-  a;u)Sjic(a  -  2u>;a  ■-  lo). 

(32) 


We  now  consider  some  cases  of  equation  (32). 


Setting  i  =  k  and  j  ~  I  and  using  the  symmetry  properties  given  by  equa¬ 
tions  (26)  and  (27),  v't  obtain  the  following  expression  for  the  asymptotic 
variaitce  of  the  cyclic  cross-periodogram, 

var{liP{a\Lj)}  ^  5„(U;  to)5jj(0;  o  -  a;)  +  |5ij(2a’  -  o;co)|‘\ 


Setting  i  =  j  and  k  ~  I  yields  the  covariance  between  two  cyclic  auto- 
period  vttgrams. 


co7;{4^\a;  w),  7^  ■'(a;  u.')}  ~  5,j(0; a;)5,j(0;  a  -  u)  +  5,j(2w  -  a;  ac)5'*,(2u7  -  a;  uj). 


liT) 


(34) 


With  i  -  j  =  k  =  I  w  e  obtain  the  variance  of  a  cyclic  auto-periodogram, 

var{lP\a\Lu)}  ^  S.,{0]ic)S„{li]Cx  -  lo)  +  |5„(2lj  -  q;u.-)1*.  ^3^) 
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suninKiri/c  these  results  in  the  asymptotic  variancc-covtuiance  matrix 
as  folio w^. 


t 

/*- 

jj 

...el 

.•.■,,(0, -l.c  .(0-1,  ..,1 

U.) 

-1-  o;.*-) 

,2 

+  S,!--  -  Oi'.'l-S'*,  (5,-  -  o,..-! 

.  .9.,(2..,  a,;...)s,,(2w 

.  .<,,(2...  -  o;,.,lSj,(2.., 

-  a.  -e) 

uO.S,(0;o 
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5.  Smoothed  L'yclic  Cross-Periodogram 


It  is  well  known,  for  the  case  of  stationary  processes,  that  the  periodograrn 
leads  to  inconsistent  estimates.  Consistent  estimates  ate  obtained  by 
convolving  tl.e  periodograrn  with  a  spectral  window,  which  amounts  to 
performing  local  averaging.  The  choice  of  window  allows  the  usei  to 
perfonn  a  bias/variance  tradeoff  in  the  resulting  estimate.  In  our  case,  as 
we  have  seen  inTheorem  1.  the  cyclic  cross-periodogram  u’)  is  an 

inconsistent  estimator  of  the  cvclic  cross- soectmm  ,9.  fo:  uj).  As  showm  bv 
Dandawave  and  Giar.nalas  [1,2],  smoothing  of  the  cyclic  peiiodogram 
results  in  consistent  estimates  of  the  cyclic  (auto)-spectrum.  Here  we 
extend  the  idea  to  the  cyclic  cross -spectrum  with  the  same  results. 

Consider  the  smoothed  cyclic  cross-periodogram  estimate  of  the  cyclic 
cross- spectrum  given  by 


c(T) 


7-1 


5=0 


i^TT^  /7n.  O.TiS. 


137) 


The  spectral  w'indow  obeys  certain  well  known  properties,  e.g., 

see  Brilliiiger  [3,  chapt.  5].  These  are  briefly  outlined  next  in  Assumption  3. 

Assumption  3 

\V{u;)  is  an  even,  real- valued  function  satisfying 


/OO  rOO 

-oc  J~oo 


(:t8) 


Tire  ^Kxtral  window  weighting  function  VV’Ti(;j)  is  constructed  from  Wlu;)  as 


where  B'pis  a  sequence  for  T=  1,  2, .  .  .  ,  with  B'f>0,  B-p--^  0,  and 
BqT  — ►  oo  as  T  ^  oo. 

We  now  present  a  theorem  that  describes  the  asymptotic  Ijehavior  of 
S’P  (a;  ui) ,  given  by  equation  (37),  demonstrating  its  consistency  and 
asymptotic  normality. 

Thaorem  2 

Suppose  X(t)  is  cyclostatic  nary  with  zero  mean  as  defined  in  section  2. 
Linder  Assumptions  1  through  3,  then 

]\m  sl]\a\Lj)  =  Sifa^u),  , 


and 

BTTam{ siP  (a;  w) ,  (a;  n) }  =  €^cav{pp  (a;  u) ,  pP  (a;  fi) } , 

where  £w  =  /.„  \W {t)'^ dr  is  the  window  energy  and  the  cyclic  cross- 
periodogram  covariance  is  given  in  Theorem  I.  Also,  the  resulting  cyclic 
spectrum  estimate  is  asymptotically  normally  distributed. 


6.  Conclusions 

Theorem  2  is  the  main  result  of  this  report  and  is  useful  for  a  number  of 
reasons,  llie  finite  b'ouricr  transfonn  X^^Xuj),  as  defined  by  equation  (5), 
c;ui  be  evaluated  using  tlic  FIT  algorithm  for  composite  record  length  T, 
so  pp(a\u)  can  be  readily  obtained.  The  smoothing  step  necessary  to 
obtain  Slj\a\ u)  via  equi.tion  (37)  is  a  sfaightforward  local  averaging 
procfxlure.  Thus,  the  smoothed  cyclic  cross-periodogram  estimate  of  the 
cyclic  |■•'•o.ss-spcctrum  can  be  readily  computed. 

The  asymptotic  variance  expression  ot  theorem  1 ,  combined  with  the 
asymptotic  normality  of  the  estimate,  allows  for  asymptotically  optimal 
criteria  to  be  developed  for  a  variety  of  detection  and  estimation  schemes. 
I'or  example,  one  can  develop  optimal  detection  tests  for  the  presence  ol 
cyclostati(.)iiary  Irehavior  and  estimate  the  cycle  frequencies  [9].  Finally, 
we  note  th:it  tlie  cioss-specUajm  and  cros.rconelation  play  a  primary  role 
in  multi-sensor  detection  and  estimation  problems.  For  example,  the  time 
dilference  ol  atrival  lietwecn  seirsors  can  be  used  u.  estimate  the  angle  of 
arrival  of  a  signal,  a  concept  that  easily  extends  to  cyclostationary  signals. 
Thus,  the  results  of  this  re|Mjrt  will  prove  useful  in  analysis  of  .such  multi¬ 
sensor  cyclostationary  signal  priKcssing  scenarios. 
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Appendix. — Proofs 


Proof  of  Lemina  1 

Using  equation  (12)  in  the  left-hand  side  of  equation  (18)  yields 


(7-1) 

Sa)  .  .  .  ,Uk-\)  - 


(i;  ,  .  .  .  ,  LJfc-l)  -  ■■  ^  <^01,  -  ,ak(i;  U- ■  ■  •  , 

ti=-(7'-i)  T*_,=-(r-i) 

y  |Cai  7")  .  •  ■  •  ,  T*._  1  )| 


rjH - -  1  7-k  I ) 


|Ti|>r  irfc_ii>7' 


E  ■  ■  E  .  -  0(T-'). 


|ti|>7-  |T;,_,i>7- 


where  the  equality  in  the  last  line  follows  after  applying  Assumption  la.D 

Proof  of  Corollary  1 

Corollary  1  follows  ininied'ately  from  Lemma  1  with  k  -2,  =  i,  and 

Proof  of  Lemma  2 

Using  equation  (5)  and  the  muliilincarity  properly  of  cumul'mls  [3,  eq 
(4.7)]  we  have 


£i=o  (*--,=0 

Using  the  assignments  t  =  4,  f]  =  4  -  '^1  =  h,~  ^  4-i  ~ 

and  A  XZ  .  uh  ,  wc  can  rewrite  equation  (43)  as 


ewri{ ^  (wi )  r  •  ■ ,  Vfc ) } 


.  (44) 


where  we  have  made  use  of  equations  (8)  and  (10).  Now,  letting  4  = 
-miti(T],  .  .  .  ,  0),  and  4  =  T--  1  -  max(T,, ....  0),  the.i  wc  may 

write, 

cixm{  X  (p  (cJi  ),•■•, 

^  E  E  Ee  (45) 

Tl-  (7  -1)  Tl,  l--(7'  l)t  £0 
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for  0  /  1  If  this  last  inequality  is  nnl  salistird  then  the  iimils 

of  the  summation  arc  taken  to  he  the  empts  set 

It  may  be  shown  that  |31 


u 


1=1. 


/  1 


t-^O 


f{. 


(16! 


where  It,  <  . Itjc- 1 1] .  with  '■  a  constant.  Using  this 

and  the  Fourier  relation 


(^  >  '  I )  *  *  *  )  "Tic  —  1  )  "y  "  ,■  -.ajt  M  j  '  '  ’  »  ^  fc  —  1  j 

a 

allows  us  to  write. 


(47) 


7’-l 

=  Ee-"' 

7-1  r-1 

E  E 

(=0 

n  =  -(T-l)  r*.i=-(7’-l) 

7-1 

7-1 

E 

E  E^«i.' 

- h^k-f'k  l) 


Cai,  --,a>;  I  7"] )  ■  ■  •  )  7fc_l) 


(48) 


ri  =  -(r-l)  t*.i=-(7’-1)  « 


The  second  term  on  the  right  hand  side  of  eciuaiion  (48)  is  bounded,  as  we 
show  next.  First,  use  equation  (47)  with  i  =  0  in  the  second  term  of  equa¬ 
tion  (48),  and  then  note  that  the  result  is  less  than  or  equal  to 


'r-i  r-i 

L  •••  Z  [|Tl|  +  -'-  +  |Tjt_l|l|Ca,,...,a*(0;Ti,---,Tfc_])| 
n=-(7'-l)  n^:=-(T-l) 

by  a  iplying  Assumption  la.  Thus  the  second  term  of  equation  (48)  is 
bounded  as  follows: 

U  ■■■  ^Cai.....a*(a;ri,  •  ,Tfc_i) 

ri  =  ^-(7’-  I)  T-fc_i=-— (T-1)  “ 

T-i  7'-l 

<  E  •••  L  [|ri|+-  -  +  |r;c-l||Cai,...,a,(0;Ti,-  -,T;fc_l)  =  C>(l), 
ti  =  -(7--1)  rfc_i=:.-(7-.l) 

and  we  can  rewrite  equation  (48)  as 
cum{A:^p(wi),  •  •  • ,  A'(p(u;*;)} 


(49) 


(50) 


ai 

r-1  7-1 

-jXt 


7'- 1 


E 


E  (f;ri, ,T,_0  f  C(l).  (.51) 


Ti=-(r-i)  Tt_.i=-(T-i) 


16 


IIk'  result  follows  by  usinji  Ixrrtma  1  in  equation  (51).L1 

Pro4jf  of  (Jor-ollary  2 

This  result  follow  s  immediately  from  I>cmma  2  with  k  -  2,  a,  =  i,  and 
«2  =  J-n 

Proof  of  Theorem  1 

First  we  show  that  uj)  is  an  asymptotically  unbiased  estimator  of 

5',j(q;  Taking  the  expectation  of  the  cyclic  cross-periodograrn  yields 

£|/;;''{a;a;)l  =  ,52) 

Next  we  employ  substitutions  similar  to  those  used  in  the  proof  of  Lemma 
2.  Let  r  =  s  minfr,  0),  and  4  =  T-  1  -  niaxfr,  0),  leading  to 

''  t-ta  T=-(r-l)  '■  ' 

Now,  using  Corollary  1,  we  obUiin, 


=  i  f;  e’'"' w)  +  0(7' 


•V 


^  (=0 

That  ill  asymptotically  unbiased  is  confirmed  by  taking  the 

limit  of  equation  (54)  as  T— >  oo,  which  results  in  Sy(a;  u). 

Next  we  develop  an  expression  for  the  asymptotic  covariance  of 
I^^^  (a;  w) .  A  theorem  of  Brillinger  [2,  theorem  2  3.2] ,  based  on  a  previ¬ 
ous  result  of  Leonov  and  Shiryaev,  provides  a  general  method  for  express¬ 
ing  a  joint  cumulant  as  a  sum  of  cumulants.  Using  this,  we  can  express  the 
covariance  of  the  cyclic  cross  periodogram  as, 

cov { lip  (a;  uj) ,  4J  ^  (/?;  /r) }  =  cum{ (u))Xj^  (a -uj),  xf  (- p)  ^  (/i  -  /3) } 

-  ~cum{ Xp (w) ,  Xf  Ho:  -  c^) ,  XlP (- /x) ,  XP'> iti-0)} 

+  ~cam{xP\u),X^k\-!^)]  cum{XP[a  -  cj),xP\^  -  /?)} 
-h  ^cum{xP(u),xPP/x  ~  p)}  cum{xP\a  -  a;),xJ^(-M)}, 

(55) 
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where  we  have  made  use  of  the  fact  that  coij{x,  y]  -  cum|x,  w*)  [2], 
Using  Lemma  2  and  Corollary  2  in  equation  (55)  yields 


=  ^  E  uj,  a  - 

''  t=o 


+ 

+ 


L  •' 


L-*  «=o 


r-i 


T 


£=Q 


T-1 


''  £=0 


+  0{T-'). 


(56) 


By  assumption,  u;,,  u>2,  cjj)  exists  and  is  finite.  Therefore,  the  first 

term  of  equation  (56)  goes  to  zero  as  co,  and  we  are  left  with 

lirn  cov{liP(a;  u),7iP(/3;  fi}}  =  S,k{uj  -  y:  uj)Sji{a  -  u  +  y  ~  (3\a  -  uj) 

^  5,<(a;  + /i  -  (5;  a;)5jfc(Q  -  w  -  u;  a  -  u;).  □ 


Proof  of  Theorem  2 

Tlie  proof  of  Theorem  2  is  somewhat  lengthy,  but  it  follows  along  the  lines 
of  Dandawate’s  Theorems  2  3^  and  2.3.5  For  this  reason  we  omit  the 
detailed  proof  here.  The  proof  of  equations  (40)  and  (41)  follows  by 
generalizing  Dardawate’s  'nieorcm  2.3.3  to  the  cross-pentxiogram  case. 
Proof  of  the  asymptotic  normality  can  be  achieved  by  showing  that  the 
higher  order  (k  >  2)  cumulants  of  the  estimate  vanish  asymptotically,  as  in 
Dandawate’s  Theorem  2.3.5. □ 
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